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in 

The first two terms in the large N asymptotic expansion of the /3 moment of the char- 
(— i ■ acteristic polynomial for the Gaussian and Laguerre /3-ensembles are calculated. This 

Q_|. is used to compute the asymptotic expansion of the spectral density in these ensembles, 

^ fH in the exponentially small region outside the leading support, up to terms o(l) . The 

leading form of the right tail of the distribution of the largest eigenvalue is given by the 
density in this regime. It is demonstrated that there is a scaling from this, to the right 
tail asymptotics for the distribution of the largest eigenvalue at the soft edge. 



> 

: 1 Introduction 

in 

One of the best known results in random matrix theory is the Wigner semi-circle law (see 
e.g. [2Ql HO])- In the case of the Gaussian orthogonal ensemble (GOE^), which consists of 
N x N real symmetric matrices H = (X + X T )/2 with the elements of X independent standard 
Gaussians, this law corresponds to the limit theorem for the eigenvalue density 

» : ssl yi'uri^o - { | (1 " A ' 2)1/2 ' \x\ > !. (") 

Lesser known is the asymptotic form of p^ Ejv (v^2iVX) in the region \X\ > 1, where its 
scaled limit vanishes. Relatively recently [15] (see also [TQl Section 1.5]) this problem arose in 
the study of the expected number A^ of critical points for the random energy 

N 

U = - ^ X j + ^( x i> -~i x n) 

where \i > and V is Gaussian distributed with zero mean and covariance 

(V(x 1 )V(x 2 ))=iV/(^(x 1 -x 2 ) 2 ). (1.2) 
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The function / in (ll.2p is arbitrary apart from the requirement that J := y/ f" (0) is well 
defined, and in fact M only depends on / through J. Thus it was shown that Af is equal to 
the product of terms smooth in fi/J times 

e -^-WJ)2)/4 p GOE JV+1( ^ 7 2 ((/i/J) + t)) ^ 

This integral exhibits two distinct behaviours depending on whether /i/J < 1 or ///J > 1. 
While knowledge of (11. ip is sufficient to determine the former, to analyze the latter requires 
the N — > oo asymptotics in the region s > 1. The density is then exponentially small, with its 
explicit form derived in [15] as 

pf° EN+1 (y/2Ns) ~ exp ( - iV^v^ 3 ! + log(s - v 7 ^!)) + 0(1)) . (1.3) 

The result (II .3p immediately raises two questions: that of deducing the explicit form of the 
0(1) term, and that of deriving the analogous asymptotics for other random matrix ensembles. 
It is the aim of the present paper to answer these questions. Motivation for pursuing this 
line of research comes from the relevance of ( 11. 3p to the right tail large deviation form of the 
probability density function (PDF) p% OE (s) for the largest eigenvalue in GOEat. 

Thus, according to [TQl eq. (14.136)] we have that for s > 1 

p%°*(V2Ns) ~ P ^ E (V2Ns). (1.4) 

To derive (11.4j) we first note that 

^ OE ( S ) = -^ OE (0;( S ,oo)) (1.5) 

where E^ OE (0; (s, oo)) is the gap probability of no eigenvalues in the interval (s, oo) of the 
GOE. But pm eq. (8.73)] the gap probability can be expressed as a series in terms of the 
fc-point correlation functions k = 1, 2, . . . according to 

pf^ E (V2Nt)dt 

(\/2N) 2 f°° f°° , / 

+ ~^- L J dh J dt 2 pg° E (v / 2iVt 1 , V2Nt 2 ) + ... (1.6) 

Assuming now that s > 1 and noting that by asymptotic independence we must have, for 
ti, t 2 , . . . , t k distinct, 

k 

pgf (V2iVti, . . • , V2Nt k ) ~ JJpgf E (v / 2iVtO, (1-7) 

1=1 

the asymptotic form (I1.3P tells us that successive terms in (I1.6P are exponentially smaller 
that previous terms. Consequently, after differentiating as required by 1 11. 5ft . ( II. 4p follows. 
Substituting ( 11. 3D in ( 11. 4p gives the leading large form of p% OE (\/2Ns) for s > 1. The latter 
was first obtained in [19] via a Coulomb gas analysis. As an application of our determination 
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of the analogue of fll.3p . and its extension by the explicit determination of the 0(1) term, for 
other random matrix ensembles we can use the analogue of (11.4j) to deduce a corresponding 
large deviations asymptotic form of the right tail of the distribution of the largest eigenvalue 
in the ensemble. 

The ensembles to be considered are the Gaussian and Laguerre /3-ensembles. They can be 
constructed out of certain tridiagonal and bidiagonal random matrices with independently dis- 
tributed elements [8]. These matrices in turn, for the special /3 values 1, 2 and 4, are Householder 
similarity reductions of the Gaussian matrices used to construct the classical Gaussian and La- 
guerre orthogonal, unitary and symplectic ensembles. For present purposes the /3-ensembles 
can be defined by their eigenvalue PDF. Generally, let us denote by MEp^w) the eigenvalue 
PDF proportional to 

N 

HwiXt) n \\k-\jf. (i.8) 

1=1 l<j<k<N 

Then the Gaussian and Laguerre /3-ensembles correspond to the eigenvalue PDFs 

ME^e"^ 2 / 2 ) and ME^A^e^x)) 

respectively, where xa — 1 f° r A true and xa = otherwise. 

Our aim is to compute the leading two terms in the asymptotic expansion of the density 
for the Gaussian and Laguerre /3-ensembles, in the large deviation regime outside the leading 
order support. This will be done by expressing the density in terms of a particular moment 
of the characteristic polynomial. The latter can be interpreted in terms of the characteristic 
function of a linear statistic, for which there is a known asymptotic formula. This then leads 
to the sought asymptotic expansion of the density. 

In Section [2] we give the details of the relation between the density and a moment of the 
corresponding characteristic polynomial, and we revise how this relates to the characteristic 
function of a linear statistic. The large N form of a large class of linear statistics for PDFs 
ME^jv(to) is expected to exhibit a Gaussian form, with explicit integral expressions for the 
mean and variance. The latter are evaluated in Sections [3] and H] for the Gaussian and La- 
guerre /3-ensembles respectively, and the leading two terms in the asymptotic expansion of the 
corresponding densities outside the leading order support are given in those sections also. A 
discussion of consequences of these results with regards to the asymptotic form of the right tail 
of the distribution of the largest eigenvalue is given in Section [51 as is a possible experimental 
realization in the case (3 = 2. 

Although the Gaussian expression for the asymptotic form of the characteristic function 
is well founded, its rigorous proof in cases corresponding to the moments of characteristic 
polynomials is known only for ensembles with two soft edge boundaries, or a hard edge with 
parameter a = [3J. This covers the Gaussian /3-ensemble, but only the a = Laguerre (3- 
ensemble. On the other hand, it is known [2] that for (3 even these moments can be expressed 
as /3-dimensional integrals. In the Appendix we show that the rigorous large N expansion of 
the latter in both the Gaussian and Laguerre cases agrees with our earlier working. 
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2 The density as a moment of the characteristic polyno- 
mial 



2.1 Exact expressions 

For the matrix ensemble ME / g i 7v(w) the eigenvalue density p(i),„(x) is given in terms of a multiple 
integral by 

P(i)A x ) = -F, — r~] w ( x ) / d\ 2 w(\ 2 )... dX N w(X N ) | [ \\k-\j\ p , (2.1) 

C PA W \ J-oo J-oo l<j<k<N 

where 

/oo /-oo 
dXiw(Xt)... dX N w(X N ) Y\ \^k-Xjf. 
00 l<j<k<N 

Defining for general g the canonical average with respect to ME^^(w) by the multiple integral 

{g(Xi, . . . , Ajv))mEsjv( w ) — r-r / dX 1 w(X 1 ) . . . dX N w(X N ) 

^p,NY w \ J-oo J-oo 

x II |A fc -A,|^(A 1 ,...,A w ), (2.2) 

l<j<k<N 

we have that 

v 

T]>-A*r) (2.3) 

corresponds to the /x-th moment of the modulus of the characteristic polynomial of the matrix 
ensemble. 

It is immediate from the definitions (12. II) and (12. 2 1) that the eigenvalue density can be written 
in terms of a moment of the characteristic polynomial, 

P(W*) = ^^% M -(x)( n i x - ^rt . (2.4) 

The corresponding inter-relation in the case of the real Ginibre matrices (N x iV matrices with 
entries independent standard Gaussians) was noted by Edelman, Kostlan and Shub [7], and 
applied for purposes of computing the density and expected number of real eigenvalues. 

For the Gaussian and Laguerre /3-ensembles we want to scale the eigenvalues so that for 
iV — y oo the support of the density is a finite interval. For this purpose let M be such that 
M/N 1. Then it is well known [9] that for the Gaussian /3-ensemble with A; i— > y/2MXi, and 
the Laguerre /3-ensemble with A; i— > 4MXi, the N — > oo support is the intervals (—1,1) and 
(0, 1) respectively. In these scaled coordinates |A| > 1 will then correspond to an exponentially 
small density for large N. 

Introducing these scalings, for the Gaussian /3-ensemble we then have 

TV 

Vm P ? 1)N+1 (Vmx) = A G e- pM ^(T[\x-x l fy 2 (2.5) 
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where 

r< r„- / 9A 2 /2l 

A G = (N + 1)(2M)^)/ 2 [ 
while for the Laguerre /3-ensemble 

N 

4M„f 1) , N+1 (4M*) = AV«V»*< £ |< - ^l^,,, (2.6) 

where 

A L = (iV + l)('4M) JV/3+1+a/3//2 C '^^[ Aa/3/2e " /3A/2 ] 

The normalizations in ( I2.5p . (I2.6P are particular Selberg integrals (see [101 Ch. 4]) and as such 
have evaluations in terms of products of gamma functions. This provides us with the explicit 
expressions 

A o _ (N + l) ( )(N0+1)/2 r(l + p/2) 

- (2tt)V2 [2Mfj} r(l + (iV + l)/3/2)' [2 - 7) 

a l _ (N+ V\(2MB) Np+1+aP/2 r (! + /V 2 ) / 28 n 

-^ + i >^ M P> r(l + (N + l)P/2)r(af3/2 + l + NP/2y [Z *> 

For future purposes we note that in the case M = N + 1 Stirling's formula applied to ( 12. 7\i 
and (12.80 shows 

y4 G ~ K( N/ 3)-m 2 (N+VWe( N+1 W 2 r(l + (5/2), (2.9) 

7T 

A L ~ Kwpyfi/^W+VWeW+Wm + 0/2). (2.10) 

2.2 Asymptotic form of the averages 

Consider the linear statistic 

N 

V(x) =^log|x-A,|, (2.11) 
i=i 

where A/ are chosen with PDF (II. 6p . The distribution of this linear statistic is given by the 
canonical average 

N 

P v (t) = (8(t-^og\x-X l \)) 
Taking the Fourier transform of both sides gives the corresponding characteristic function 

P v (k) = ( e ik ^ l0 ^ x - x A . (2.12) 

Observing 

N 



i=i 



ME 0iN (w) \ I ME^jv(iii) 
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we thus have 

N 

T]>-A/\ =iVH/3). (2.13) 

More generally let G := Ylf=i9(^j) denote an arbitrary linear statistic. The characteristic 



function is then 



N 



p G (k) = ( TT e ^E i=1 9^)\ _ / 214) 

The significance of the relation (12.131) is that for a large class weights w and statistics g, 
(12.141) has a known asymptotic form. Thus suppose that for large N the eigenvalue support of 
ME/3 ]sr(g) is the finite interval [a, b]. Suppose too that g(X) and its derivative are bounded on 
[a, b] or more generally that (I2.17P below is finite. In these circumstances there are well founded 
physical grounds (see e.g. [TQ], Section 14.4.1]) to expect that Pc(k) is a Gaussian, 



Jk m (g) e -k 2 (a(g)f/2 

where the mean Hn{9) and variance (a(g)) 2 have the explicit forms 



P G (h) ~ e ik m (g) e -k^*(g)y/2^ ( 2 _ 15 ) 
N-+oo 



M9) = j P(i),N(t)g(t)dt (2.16) 

a{9)) =wL dtl w^m^wh dh ^ 

hf^kal a k :=\ f g + ^ cos o) cos k9 d9 (2.17) 

' 7 1 



k=l 



(see also [H]). At a rigorous level, for a class of weights w(x) including the Gaussian but not 
Laguerre, and a class of g(x) which does not include g(x) = log \x — t\, this has been proved 
by Johansson [16]. Very recently [3], for matrix ensembles in which the eigenvalue support is 
a single interval, with both endpoints soft edges or a hard edge with parameter a = 0, the 
conditions on g(x) have been relaxed to require only that g{x) be analytic in a neighbourhood 
of the support. The Gaussian /3-ensemble has the first of these properties, with g(x) = log \x— 1\ 
for x outside the eigenvalue support has the second of the properties, and thus then (12 . 151) is 
a rigorous theorem. But in general the Laguerre ensemble has a hard edge with a ^ 0, and so 
the result of [3] applies only in the special case a = 0. 

We will proceed under the assumption that (" 12 . 1 5[) is valid for the linear statistic (12. lip not 
only in the Gaussian case, but also in the general Laguerre case, with x outside the interval of 
support of the eigenvalue density [a, b]. Then use of f 1 2 . 1 5 j) in (I2.13P gives 




with v{t) = log \x — t\ and appropriate w. It thus remains to evaluate /iAr(f) and cr{v). Sub- 
stitution in (12.51) will then give the sought large N form of the densities. We will treat the 
Gaussian and Laguerre cases separately. 
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3 The Gaussian /3-ensemble 

For the weight w(X) = e~ M/3x2 in ( 11. 8ft with M/N — > 1 as N — > oo the leading form of the 
density p(i),jv(t) is given by the RHS of (II. ip multiplied by N. The next order term consists of 
an oscillatory and a non-oscillatory part jl] - - only in the latter contributes to the integral in 
(I2.16P to the same order. Ignoring then the oscillatory term we read off from Johansson [T6t 
displayed equation below (2.15)] or more explicitly [6] eq. (19), Lemma 2.20], [121 eq. (1.4)] 
that for \t\ < 1 

, , 2(MN) 1 ' 2 / M 2 \V2 /l 1\/1 P/ ,x I,, ^ 1 1 



7r V N J \/3 2/V2 v ' 2 v ' ir^/T~fi 

2M (1 oM/2 iN/lr/ ,s l r/ I I 

;i_^)i/ 2+ ) -5(t-i) + 



vr v ' V/3 2A2 v ' 2 v ' vrv^3t2 
N-M 

where the second asymptotic equality follows by expanding M/N about l. 

We seek the value of (I2.16P with P(i),jv(t) given by 03. H) and a(t) = log(x — t). For this we 
make use of the following integral evaluations. 

Proposition 1. Suppose x > 1. We have 

of 1 1 
- / log \x - t\ (1 - t 2 ) 1/2 dt = x 2 - xVx 2 - 1 - log(2(x - Vx 2 - 1)) - -, (3.2) 
7T 7_i 2 

t 2 )- 1 /2dt = iog(i(a:+Vi Ji: T)). (3.3) 

Proof. It is straightforward to check that for both equations, both sides tend to logx + 0(1 /x) 
as x — > oo. Hence it suffices to show that both sides have the same derivative with respect to 
x, and thus that for x > 1 

- / 11 t} rft = 2x(l-(l-l/x 2 ) 1/2 ), 
J-i x — t 











/jog 


x — t 


















\x — t 


* 







7T 



i r 1 i dt i 



7T j_ x {l-t 2 ) l / 2 X-t (X 2 - l) 1 ' 2 ' 

The first is a well known result in random matrix theory (see e.g. pTOj. eq. (1.32)]). Both 
integral evaluations can be checked by expanding the left hand and right hand sides in powers 
of 1/x, and making use of the beta integral on the LHS and the binomial expansion on the 
RHS. □ 



Making use of Proposition CD we therefore have 

fi N =m(x 2 - xVx 2 - 1 - log(2(x - Vx 2 - 1)) - - 



X 2 



+ (N-M) log (-(x + Vx 2 - 1)) . (3.4) 
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Next we are next faced with the task of evaluating (" 12. 1 7[) for a(t) = log \x — t\ and a = — 1, b = 1. 
Making use of the second equality in (12.171) . simple manipulation shows we must evaluate 

a k = - f log ( 1 - — ) cos k9 d9 (3.5) 



TC J \ X 

for k = 1, 2, . . . and x > 1. 
Lemma 2. For s > 1 

log (l - — ) = log(l - ue l6 )(l - ve-' s ) - log(l + u 2 ), (3.6) 

where 

U = X - ( X 2 - 1)V2 (3 ^ 

and thus in particular < v < 1 . 
Proof. The RHS of (I3.6P is equal to 

Equating with the LHS we read off that 

1 2u 



x 1 + v 2 

and solving for v gives (13.71) . □ 
Substituting (13. 6 \ in (13. 5 \ shows 

a k = — — , (3.8) 
and substituting this in the second equality of (12.171) we see upon evaluating the sum that 

o- 2 = -|log(l-^) 

= -- log(a; 2 - If I 2 - | Iog(2(a: - Vx^~l)), (3.9) 

where the second equality follows upon making use of (13. 7\i . 

Let us now make the choice M = N + 1 . Then it follows from (I2.18P , (13. 4 p and (I3.9P that 

N \ 

(N+l)f3(x 2 -xVx T ^l-\og{2{x-Vx I: ^))-l/2) 



\\\x-\i 



/ W-5-oo 
>*=1 / ME a , JV (e-( JV + 1 )' 9A2 ) 



( x 2 _ 

(3.10) 



((X+ (X 2 - l) 1 /2)/2)(l-/3/2)' 



We remark that with the exponent in the average on the LHS replaced by 2a, and /3 in ME^jv 
set equal to 2, the asymptotic formula (12.151) has been used in jTTj to deduce the corresponding 
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1.2 


1.4 


1.6 


1.8 


1.2 


1.4 


1.6 


1.8 


6 


1.404 


1.177 


1.109 


1.078 


1.126 


1.074 


1.056 


1.047 


12 


1.220 


1.092 


1.055 


1.039 


1.067 


1.038 


1.028 


1.023 


18 


1.152 


1.062 


1.037 


1.026 


1.046 


1.025 


1.019 


1.015 


24 


1.116 


1.047 


1.028 


1.019 


1.035 


1.019 


1.014 


1.011 


30 


1.094 


1.038 


1.022 


1.015 


1.028 


1.015 


1.011 


1.009 



Table 1: The ratios ^pjy(s) (first four inner columns) and r^ N (s) (final four columns) for varying 
iV (rows) and given s (columns). 



asymptotic form up to the same order as in (13.101) . With a = 1 the result of [T7] agrees with 
the = 2 case of fl3TT0|) . 

The result (13.101) substituted in ( 12 .51) and use of (12. 91) then gives our sought generalization 
of ([I3D, 



2Npf l)N (V2Nx) ~ e-A^^T+iog^-^I)) 

iV— >oo 



valid for x > 1. As a check on our working we note that in the case /3 = 2, corresponding to 
the GUE, the leading two terms in the large N expansion of the largest eigenvalue distribution 
p^ UE (v2iVs), s > 1, has been obtained recently [2TJ eq. (76) with a = 1 and i = \/2iVs]. But 
the analogue of (11.4j) for the GUE tells us that this asymptotic form must be identical to that 
for pSx N (\/2Ns) 1/3=2- Indeed ( 13. lip with /3 — 2 reproduces the result from [2T] . 

In the cases (3 = 1,2 and 4 simple closed form expressions for p£L ^(x) are known in terms 
of Hermite polynomials (see e.g. [101 Ch. 5&6]). The simplest case is (3 = 2 for which 

pfi)A x ) = ^ N _ 1 y\ H N(x)H N -i(x)-H' N _ l (x)H N (x)y (3.12) 
while for /3 = 1 and even we have 

-x2/2 X 



P(i),n( x ) 



e~ ' 1 /"°° 2 



(3.13) 

The case /3 = 4 is similar to (I3.13p . but to make our point is suffices to restrict attention to 
P = 1 and 2. 

Numerical evaluation of (13. 12[) and (13 . 1 3[) allows the accuracy of (13. lip to be quantified. 
Define the ratio 



tAs) (3,4, 



where p^ a ^ m (v / 2iVs) is the asymptotic form ( 13. lip . The results for particular s and varying 
are given in Table [U 
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4 The Laguerre /3-ensemble 

We now turn our attention to the weight w(X) = A a/3//2 e~ 2M ^ A XA>o in (I1-8D with M/N — > 1 as 
iV — > oo. The leading two terms in the large N expansion of the smoothed density for < t < 1 
is most naturally presented in the variable t 2 . Thus according to [T2l eqns. (6.21) and (6.22)] 
we have 

, 4(MA^) 1 /2, M 2 \V2 / 1 In/ 1 1 , . 



1 1\ /l P/ 
B -2) (2^ 



+ V/3 2/V2'" t\/T~7? 
4M 



7T 



^M^-20(^->) 



VT 



+ -.) ' ) + « (4.1) 



Noting that 



y \og\x -t\p w<N (t)dt = J log 



x-t|(2tp ( i )>JV (0)d*, 



where we define 2tp(^ N {t 2 ) for £ < by the requirement that it be even (equivalently replace 
S(t — 1) in (14. ip by 5(t — 1) + 5(t + 1)) we see that the corresponding form of f l 2 . 1 6 [) with 
a(t) = log \x — t\ can be simplified using knowledge of (I3.14p and (13. 3p . Thus, with M = N + 1 



fJL N 



--2M(x - ^x{x - 1) - log(2( v / x - y/x - 1)) - -) 
+ (a - 2) log (-(V^ + V^ 3 !)) - ^(2a + 2//3 - 1) log |x| 



1/2 



+ " 2) log ( X - X ) 1/2 - ^ 

Our remaining task is to evaluate (12.171) for a(t) = log \x — 1\ and a = 0, 6 = 1. Making use 
of the second equality in (I2.17P shows we must evaluate 



Wo H 2x-lJ 



for k = 1,2,... and x > 1. For this we can use Lemma [2] with x replaced by 2x — 1 in (13. 7p 
and thus 

v = 2x - 1 - 2x 1/2 (x - 1) 1/2 . (4.3) 

In particular it follows that is again evaluated according to (13.81) . but with v now given by 
(14. 3p . Consequently 

a 2 = -|log(l-, 2 ) 

= -1 \og{x{x - l)) 1 / 2 - i log(2(x 1 / 2 - (x - l) 1 / 2 )). (4.4) 
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The results (14. 2[) and (14.41) substituted in (" 12 . 1 8[) give for the asymptotic form of the 
moment of the characteristic polynomial for the Laguerre /3-ensemble, 



N . 

^ e 2(JV+l)/3(x- v / :C ( :C -l)-log(2( v ^- v ^T)-l/2) 

N^oo 



( r _ l)(l-3/5/2)/2 

x 2 "° M g((>f U (V^ + V^T)^. (4.5) 



We remark that for fixed iV and a? large the LHS tends to or . It is straightforward to check that 
this is the x large behaviour of the RHS, suggesting that (14.51) holds uniformly for x > 1 (the 
corresponding Gaussian asymptotic formula (I3.10P also exhibits this property). Substituting 
this in (12. 6 p with M = N + 1 and making use too of (I2.10p we obtain for the asympotic form 
of the density 

ANp\ llN (ANx) ~ e -2^(v^T)+io g (v^-^T)) 



W->oo 



x (p/2+l)/2 



.(^ + v^T)^ ^-^r (l + f ) ^ W)" 



-/8/2 



(4.6) 



In the Laguerre case, for — 1, 2 and 4 there are simple closed form expressions for the 
density in terms of Laguerre polynomials (see e.g. pUl Ch. 5&6]. The simplest case is (3 = 2, 
for which Y(N ) <i d 

The explicit formulas for — 1 and 4 has a structure similar to (I3.13p . We have used these 
formulas (in the cases — 2 and 1) to check that that indeed for given a, a fixed x and 
increasing iV the ratio of the LHS and RHS in (14.61) tends to unity. 



5 Distribution of the largest eigenvalue 

We have noted in the case of the GOE that the large deviation form of the right tail of the PDF 
for the distribution of the largest eigenvalue is, up to exponentially small corrections, equal to 
the corresponding large deviation form of the density (recall eq. (I1.4p ). The given derivation 
of this formula applies equally as well to the right tail of the largest eigenvalue PDF for the 
Gaussian and Laguerre /3-ensembles. Thus with s > 1 we have 

p%(V2Ns) ~ pf llN (V2Ns), (5.1) 

where the RHS has the explicit form given by (13.111) . and 

p L N (4Ns) ~ pl )jN (4Ns), (5.2) 

where the RHS has the explicit form given by (14. 6p . 
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In [19] the leading asymptotic expressions 



p%(V2Ns) ~ exp (-/3iVV?(V2(s - 1))) , (5.3) 
p^(4iV S ) ~ exp (-»J(4( S - 1))) , (5.4) 

valid for s > 1 were derived. Here, with G(z) := 3 F 2 ({1, 1, 3/2}, {2, 3}; z) a particular hyper- 
geometric function, the rate functions ^+ and iflt are specified by 



tf{z - V2) := - log(,V2) + ^ G (|) 



(5.5) 



^(z _ 4 ) := - log* + - G ( - ) . (5.6) 



z \ z 



Subsequently in [211 eq. (15)] the first of these was written in the simplified form 

^{z -V2) = + log ^ ^= j . (5.7) 

Substituted in (I5.3P this implies 



p%(V2Ns) ~ e-^-v-'-ij+tagi-v-'-iH (5. 8 ) 

which indeed is the leading asymptotic behaviour as predicted by (15. ip and (13. lip . It is remarked 
in [5T] that (15. 8 p in the GOE case /3 = 1 was first established in an earlier work [Tj. 

Note that comparison of (15.51) and (15. 7p implies an explicit form for G(2/z 2 ). This substi- 
tuted in (I5.6P implies 

<$Xz - 4) = ^ z ( z -4) + 2 log V 7 ^) • (5.9) 

Now substituting (15.91) in (15. 4p we obtain the explicit large derivation formula 



p h N (4Ns) ~ e -2^( v / ^)+i° g (v^-vCT))_ (5 1Q) 

And as in the Gaussian case, this indeed agrees with the leading term as implied by our 
corresponding formulas (15. 2p and (I4.6p . 

The large deviation tails of the PDF for the distribution of the largest eigenvalue in the (3 = 
2, a = Laguerre ensemble have been the subject of a recent experimental study [2]. We note 
that this particular Laguerre ensemble gives the eigenvalue PDF for the matrix product X'X, 
where X is an NxN standard complex Gaussian matrix (see [TUl §3.2]). The experimental study 
relates to the probability distribution of the combined output power from a certain coupled array 
of high gain lasers. It turns out that the round trip propagation matrix M determining the 
evolution of the complex electric field in this setting can be decomposed in the form X^X for X 
a square random matrix with near Gaussian complex entries. Moreover, the coupled lasers will 
lase at the mode with minimal losses, which corresponds to the eigenvector of M with largest 
eigenvalue, and the output power is proportional to this eigenvalue. Significantly, very large 



12 



data sets can be generated (order half a million measurements) allowing for the large deviation 
regime of the PDF to be probed. 

In regards to the right tail, the experimental data was compared against the leading form 
(|5.10p . modified by multiplication by a numerical factor C2 = 0.0063 chosen on the basis of the 
modification then giving a better fit of results obtained from numerical simulation. But our 
result (15. ip extends (I5.10p to give the algebraic and constant terms. Explicitly, for (3 = 2, a = 
we have 



p L N (ANs) ~ e-mV^+iosiVs-V^D) 1 ( _J_ ] . (5.H) 
yNK J s s - 1 V 32vr iV/ v ; 



It thus remains to see if the right tail of the experimentally determined PDF can be resolved 
to the extent that this correction term can be exhibited. 

The formulas (15. ip and (15. 2p are valid for N — > oo with s > 1 fixed. This is to be contrasted 
to the soft edge scaling limit in which N — > oo and simultaneously s — > 1 according to the 
scalings [9] 

V2Ns = V2N + _ , , (5.12) 

ANs = AN + 2(2N) 1/3 X (5.13) 

in the Gaussian and Laguerre cases respectively. For (3 even the explicit soft edge scaling form of 
the density p s ,°^(X) — which as a consequence of the principle of universality [18] is independent 
of the particular case, Gaussian or Laguerre — has been calculated in [U Corollary 12]. Of 
interest to us is the corresponding right tail asymptotics deduced from this exact expression [U 
Corollary 14] (corrected by the removal of a spurious factor of 1/2 on the RHS) 

P W (X) xZoo n (W/2 IW4-i/ 2 ' ( 5 - 14 ) 

In the context of largest eigenvalue distributions, the argument leading to (ll.4p tells us that 
the PDF, p soft (X) say, of the soft edge scaled largest eigenvalue of the Gaussian and Laguerre 
/3-ensembles must have the same leading large X asymptotic form as p s ^(X). Thus for even 
(at least) 

i r(i + g/g) e- 2 ^ 3/2 /3 

The cases (3 = 1 and 2 of this expression, deduced using the exact Painleve transcendent 
evaluations of p soit (X) in these cases (see e.g. [TQl Ch. 9]) were reported in [22]. The equivalent 
asymptotic expression 



X 



1 r(l _L g/O) p -2/3X 3 / 2 /3 



up to the explicit form of the constant but applying for general (3 > 0, has recently been derived 
in [5], using the stochastic differential equation characterization. 
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It is possible to use our explicit large deviation forms A3. 11 j) and (14. 6[) for p£L N (y/2Ns) and 
p\i) n(^Ns) to reclaim (I5.14p . Thus as done in [19] for the leading asymptotic forms (15.81) and 
(15.101) . we replace s by (I5.12p (Gaussian case) and (I5.13P (Laguerre case). Taking the limit 
iV — > oo then reclaims the RHS of (15.141) . In fact such correction formulas to the X — > oo form 
of p^*(X) should hold for the complete asymptotic series of p9^ N (v2Ns), pV^ N (ANs) upon 
this double scaling. 

The results of our study also have consequences with regard to the large deviation right 
tail form of the probability E^^(k; (s, oo)) that there are precisely k eigenvalues in the interval 
(s, oo) of the Gaussian and Laguerre /3-ensembles. Now analogous to (II. 6p 



E N> p(k; (s, oo)) = J dti---j dt k p(k)(ti, ■ ■ ■ ,t k ) 

fk + l\ f°° , f°° , , . 

— I ^ J J dt 1 ---J dt k+ i P(fc+i)(ti, . . . , t k+ i) + . . . 

For the Gaussian case, changing variables as in (II. 6p and making use of the analogue of (JTTTJ) 
shows that for s > 1 

/ roc \ k 

fu. U /o AT ^ _ , I , /o l\T\k / / J3 



E^(k; (V2Ns, oo)) ~ (V2N) k p^ N (V2Nt) dt 

with y/ 2Np9 l -, N (\/2Nt) in turn given by (13. lip . Similarly, in the Laguerre case, we have that 
for s > 1 

(roo \ k 

j pf llN (ANt)dtj 

with ANpf 1)iN (ANt) given by (14. 6p . The soft edge scaling limit of either of these reads 

Ef\k;(s, oo))^ p s $(t)dt) . (5.17) 

J S 

This asymptotic form compliments that for E s ^ ft (k; (s, oo)) with s — > — oo recently derived in 
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A Appendix 

Here the results ( 13 . 1 Of) and (14. 5 p for (3 even will be derived by applying the steepest descent/ 
stationary phase method of asymptotic analysis to duality formulas giving the averages as 
0- dimensional integrals. 
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A.l Gaussian case 

According to PH (13.162)] 

G Vj A = (f[(c-iy) N \ • (A.l) 

' ME 2/QiJV (e-« 2 ) \i =1 ' ME 2Q , n (e-« 2 ) 

Let (3 be even and write 1/a = (3/2, n = (3. After scaling the variables (lA.lj) then reads 
In words ( 1A.2I) tells us that the (3 moment of the characteristic polynomial in ME^ 

^(e-W) i s 

proportional to the N-th moment of the characteristic polynomial at — ix in ME 4 / / 3 ( g(e~ 2My ). 
The latter is well suited to an N — > oo asymptotic analysis. 
As a multidimensional integral the RHS of ( 1A.2j) reads 



N 
\3=1 



1 



2 CV/j^[e-W 

ME 4/ft(i (e-«» 2 ) 



X 



/oo /-oo _ 

^- ^ e - 2 ^-^^=^l e ^Ef =1 /(^) JJ \ yj - yk \W (A.3) 



l<<k</3 



where 



f(x,y) := -2y 2 + log(x-iy). (A.4) 



For N — y oo this can be analyzed using the method of stationary phase 
To apply this method, we note 

|f = - — , (A.5) 
ay x — iy 

dy 2 (x — iy) 2 ^ 

The first of these implies that for x > 1 the stationary points of / occur at 

I/± = -I(x±(x 2 -l) 1 /=»). (A.7) 

We would like to deform the contours of integration from the real line as in ( 1A.3I) to pass through 
one of these points. The details of how to rewrite ( 1A.3I) so the integrand is an analytic function 
is given in [H Lemma 1]. Noting from f ]A.6j) that d 2 f/dy 2 \ _ < suggests we translate 
each contour parallel to the real axis to pass through y_ (the other choice y = y + is such that 
d 2 f/dy 2 > 0, which would mean deforming the contour to run parallel to the imaginary axis 
in the complex plane, causing convergence problems at infinity). Furthermore, we expand 



f&, Vj) = f(x, yJ) - ^a x {y- - Vjf + 
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where 



f{z, y.)=x 2 - x(x 2 - if' 2 - i + log { ^ + {X \ 1)V2 ) , (A.8) 
a x = 8(x 2 - l) 1/2 (x - (x 2 - 1) 1/2 ). (A.9) 



Substituting in (1A.3I) and changing variables shows, after setting M = N + 1 



yj=1 1 ME (3i4//3 (e- 2 ( w + 1 )w 2 ) 



Substituting flA~7|) - f[Q|) reclaims f lXTUj) 



A. 2 Laguerre case 

The analogue of (1A.2|) in the Laguerre case is the identity [101 eq. (13.44), Exercises 13.1, q.4] 



C N ,p [x^e-W] (AM) Nm I" 

\J=1 / MEjV(3 ( :r/ 3a/2 e -2Af/3 I ) 



1 I dz ± I dz m TT -AMt Zl -JV-l-(2//3)(m-l) n Na+iV 

M m (a,N,2/(3)f2ri-"f 2m ^ 



C<fc<m 

where d := a — 1 + 2//3 and 



x 

K<fc<m 



m— 1 



A/f ( R \ TT r(l + a + /3 + j 7 )r(l + (j + l)7) , A 9 ^ 

M m a,/3,7 =11— — , — — - — r. A.2 

f = o r(i + a + j7)r(i + p + n)T{i + 7) 

With m = (3 (P even) this expresses the moment as a /3-dimensional integral, and furthermore 
the latter is well suited to asymptotic analysis. But before doing so, we apply Stirling's formula 
to the known gamma function evaluation of the normalization on the LHS of (1A.1[) as used in 



the derivation of (12.101) to obtain their asymptotic form. This allows us to replace the LHS of 
(OCT]) by 



e^2 2 ^r(a/j/2 + l)r((a + l)/3/2) « _ \ 

(nN(3)(NP/2f a+1 ^ 2 \H V „, ' 1 J 

v yj \ y/ j \j=i i ME w (^»/2 e -2M^) 



We also simplify the prefactor on the RHS. Thus it follows from ( IA.2j) . the duplication formula 
for the gamma function 

,3/2-1 

Yl F(z + 2j/p) = (2tt) W2 - 1)/2 (P/2) 1/2 -^ /2 T(Pz/2) (A.4) 

j=0 
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and Stirling's formula that for iV — > oo 

, 1 , , ~ (T(l + 2/ ft? mfm/2f^ r W 2 + l)rQg(a + l)/2 + 1) 
Mp(a,N,2/p) Ky ,HJJ KHI 1 KH 11 r(/9/2 + l)r(/3 + l) V ; 

We now turn our attention to the integral on the RHS of (lA.lj) . and for convenience specialize 
to the case M — N + 1, although our final formula will be valid for any M — N fixed. This 
integral can be written 

f S n eM9(Mi)(i + ^-^v™^ n & - ^) 4//3 ' ( a - 6 ) 

i=i i<i<fc</3 

where 

#(t, 2) := -4te - log 2 + log(l + z). (A.7) 
In preparation for applying the method of stationary phase, we note that 

d 9 , 1 1 



<9z z 1 + z 

3 2 g _ 1 1 

a? ~ ^2 ~ (i + 2)2- 

Thus there are stationary points at 

z± = J(-i±(i-iA) 1/2 ) 



and at these points 



7i := (4t) 2 (l - 1/t) 1 / 2 . (A.8) 



<9z 2 

Choosing the positive sign allows the contours (which must all encircle the origin) to be de- 
formed to pass through z + parallel to the imaginary axis; the choice z_ would require deforming 
the contour along the (negative) real axis, which is a branch cut of the logarithm function. Ex- 
panding about z + we have to second order 

1 



g(t, h) = g(t, z+) - -7* xj , Zj - z + 



2 i" ~Ji "J r --"J) 

with 

g(t, z+) = 2t(l - (1 - 1/t) 1 / 2 ) - 2 \og(Vt - y/t=T) 

and 7t as in (IA.8|) . Substituting in ( ]A.6j) and making a further change of variables yP%Xj = yj 
shows that for A" — > 00 (1A.6[) has the asymptotic form 

( 2 ^)V /3 , 4//3 [e^ 2 / 2 ](l + , + )^- 2 + 2 ^)^ 2 ^ 1 )(iV 7t )(2-3/3)/2 

= ^^g^^^ Wg) 1 -^^^ (1 + ^^z^ 1 ^)™'*, (A.9) 

where the equality follows upon making use of the explicit gamma function evaluation of 
Cp A /p[e~ x2/2 } [El eq. (4.140)] and the duplication formula ( ]A.4j) . 

Multiplying ( 1A.9|) and ( 1A.5I) . equating to ( 1A.3I) and simplifying reclaims ( 14. 5ft . 



17 



References 

[1] G. Ben Arous, A. Dembo, and A. Guionnet, Aging in spherical spin glasses, Prob. Th. 
Rel. Fields 120 (2001), 1-67. 

[2] T.H. Baker and P.J. Forrester, The Calogero- Sutherland model and generalized classical 
polynomials, Commun. Math. Phys. 188 (1997), 175-216. 

[3] G. Borot and A. Guionnet, Asymptotic expansion of beta matrix models in the one-cut 
regime. la7XTv:1107.11671. 2011. 

[4] P. Desrosiers and P.J. Forrester, Hermite and Laguerre (3-ensembles: asymptotic correc- 
tions to the eigenvalue density, Nucl. Phys. B 743 (2006), 307-332. 

[5] L. Dumaz and B. Virag, The right tail exponent of the Tracy- Widom-beta distribution, 
larXiv:1102.4"8T8l 2011. 

[6] I. Dumitriu and A. Edelman, Global spectrum fluctuations for the (3-Hermite and (3- 
Laguerre ensembles via matrix models, J. Math. Phys. 47 (2006), 063302(2006). 

[7] A. Edelman, E. Kostlan, and M. Shub, How many eigenvalues of random matrix are real?, 
J. Amer. Math. Soc. 7 (1994), 247-267. 

[8] A. Edelman and B.D. Sutton, From random matrices to stochastic operators, J. Stat. Phys. 
127 (2006), 1121-1165. 

[9] P.J. Forrester, The spectrum edge of random matrix ensembles, Nucl. Phys. B 402 (1993), 
709-728. 

[10] , Log-gases and random matrices, Princeton University Press, Princeton, NJ, 2010. 

[11] P.J. Forrester and N.E. Frankel, Applications and generalizations of Fisher-Hartwig asymp- 
totics, J. Math. Phys. 45 (2003), 2003-2028. 

[12] P.J. Forrester, N.E. Frankel, and T.M. Garoni, Asymptotic form of the density profile for 
Gaussian and Laguerre random matrix ensembles with orthogonal and symplectic symme- 
try, J. Math. Phys. 47 (2006), 023301. 

[13] P.J. Forrester and N.S. Witte, Asymptotic forms for hard and soft edge general (3 condi- 
tional gap probabilities. larXiv:1110.4284[ 2011. 

[14] M. Fridman, R. Pugatch, M. Nixon, A. A. Friesem, and N. Davidson, Measuring maximal 
eigenvalue distribution of Wishart random matrices with coupled lasers, larXiv:1012 .1282, 
2010. 

[15] Y.V. Fyodorov, Complexity of random energy landscapes, glass transition and absolute 
value of spectral determinant of random matrices, Phys. Rev. Lett. 92 (2004), 240601, 
Erratum: Phys. Rev. Lett. 93 (2004), 149901. 



18 



[16] K. Johansson, On fluctuation of eigenvalues of random Hermitian matrices, Duke Math. 
J. 91 (1998), 151-204. 



[17] I. Kravotsky, Correlations of the characteristic polynomials in the Gaussian unitary en- 
semble or a singular Hankel determinant, Duke Math. J. 139 (2007), 581-619. 

[18] A. Kuijlaars, Universality, The Oxford Handbook of Random Matrix Theory (G. Akemann, 
J. Baik, and P. di Francesco, eds.), Oxford University Press, Oxford, 2011, pp. 103-134. 

[19] S.N. Majumdar and M. Vergassola, Extreme value statistics of eigenvalues of Gaussian 
random matrices, Phys. Rev. Lett. 102 (2009), 060601. 

[20] M.L. Mehta, Random matrices, 2nd ed., Academic Press, New York, 1991. 

[21] C. Nadal and S.N. Majumdar, A simple derivation of the Tracy- Widom distribution of the 
maximal eigenvalue of a Gaussian unitary random matrix, J. Stat. Mech. 2011 (2011), 
P04001. 

[22] M.G. Vavilov, P.W. Brouwer, V. Ambegaokar, and C.W.J. Beenakker, Universal gap fluc- 
tuations in the superconductor proximity effect, Phys. Rev. Lett. 86 (2001), 874-877. 



19 



